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Abstract. We give a detailed description of the implementation of a M0lmer- 
S0rensen gate entangling two '*''Ca+ ions using a bichromatic laser beam near-resonant 
with a quadrupole transition. By amplitude pulse shaping and compensation of 
AC-Stark shifts we achieve a fast gate operation without compromising the error 
rate. Subjecting different input states to concatenations of up to 21 individual gate 
operations reveals Bell state fidelities above 0.80. In principle, the entangling gate 
does not require ground state cooling of the ions as long as the Lamb-Dicke criterion 
is fulfilled. We present the first experimental evidence for this claim and create Bell 
states with a fidelity of 0.974(1) for ions in a thermal state of motion with a mean 
phonon number of n ~ 20(2) in the mode coupling to the ions' internal states. 



Submitted to: New J. Phys. 
1. Introduction 

Building a device that is able to carry out arbitrary calculations by exploiting the laws 
of quantum physics has been an experimental challenge for more than a decade now. A 
large variety of physical implementations have been conceived to meet the requirements 
for quantum information processing summarized in [1]. Among these implementations, 
strings of ions stored in linear Paul traps and manipulated by laser pulses have proven 
to be a particularly successful architecture to realize quantum information processing. 
Experiments with trapped ions have shown long relevant coherence times [21 El H], 
the ability to faithfully initialize and read-out qubits O [6] and high fidelity quantum 
operations [3, [U [9]. Current efforts are focussed on scaling up ion trap experiments 
to handle many ions, improving the quality and speed of the basic operations and 
integrating the various techniques into a single system. Concerning the basic operations, 
the realization of universal multi-qubit gates is particularly challenging. Many different 
types of gates have been proposed over the last years and several of them have been 
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experimentally investigated. Gates using a collective interaction P [10] between the 
ions and the laser field - until recently only applied to qubits encoded in the hyperfine 
structure (hyperfine qubits) - were very successful in creating multi-particle entangled 
states and demonstrating simple quantum error correction techniques 

Recently we demonstrated the first application of a M0lmer-S0rensen gate operation 
to an optical qubit, i.e. a qubit encoded in a ground and a metastable state of ^"Ca"*" 
ions, deterministically creating Bell states with a so far unmatched fidelity of 0.993(1) 
[9]. Here, we present a further investigation of this universal gate operation acting on 
optical qubits and extend the theoretical and experimental analysis. Particular emphasis 
is put on the compensation of AC-Stark shifts and amplitude pulse shaping to reach high 
fidelities without compromising the gate speed substantially. The gate characterization 
is extended further by investigating the fidelity decay for different input states after up 
to 21 individual operations. 

Moreover, we report on the first experiments demonstrating a universal entangling 
gate operating on Doppler-cooled ions. We derive simple expressions [12] for the qubit 
populations under the action of the gate and and use these equations to infer the mean 
vibrational quantum number n of the axial center-of-mass mode. For ions in a thermal 
state with n = 20(2), we obtain Bell states with a fidelity of 0.974(1). 

The ability to implement high fidelity multi-qubit operations on Doppler-cooled ions 
is of practical interest in ion trap quantum information processing as the implementation 
of quantum algorithms demands several techniques that do not conserve the ions' 
vibrational quantum state: (i) State detection of ancilla qubits as required by quantum 
error correction schemes [13] can excite the ion string to a thermal motional state close 
to the Doppler limit because of the interaction with the laser inducing the ions to 
fiuoresce. (ii) Experiments with segmented traps structures where ion strings are split 
into smaller strings also tend to heat up the ions slightly [T3]. Here, the availability of 
high-fidelity gate operations even for thermal states may provide a viable alternative to 
the technically involved re-cooling techniques using a different ion species [T5l [T6] . 

2. Theoretical gate description 
2.1. M0lmer-S0rensen gate 

A two-qubit quantum gate that is equivalent to a controlled-NOT gate up to local 
operations is achieved by the action of a Hamiltonian if oc cr„ ® ct„, where (t„ = cr ■ n 
is a projection of the vector of Pauli spin matrices onto the direction n [17] . Two 
prominent examples of this type of gate are the conditional phase gate [H |T8] and the 
M0lmer-S0rensen gate [IHl EHl [10] . In the latter case, correlated spin flips between the 
states It) IT) ^ li)li) si-i^d lT)li) ^ li)IT) ^'^^ induced by a Hamiltonian 

H oc (T0®cr0 where = cos 0(7^. + sin</)crj^. (l) 

The unitary operation U = exp{ija^ ® cr^) maps product states onto maximally 
entangled states. In 1999, the proposal was made to realize an effective Hamiltonian [191 
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Figure 1. (a) M0lmer-S0rcnscn interaction scheme. A bicliromatic laser field couples 
the qubit states ^ |||) via the four interfering paths shown in the figure. Similar 
processes couple the states |ti) ^ lit)- The frequencies uj± of the laser field are 
tuned close to the red and the blue motional sidebands of the qubit transition with 
frequency loq, and satisfy the resonance condition 2ll}q = uj^ + The vibrational 
quantum number is denoted n. (b) Level scheme of '*''Ca"'"showing the transitions used 
for cooling/detecting, repumping and resetting the state of the ion as well as the qubit 
transition. The qubit is encoded in the metastable state jf) = |Z)5/2,w = 3/2) and 
the ground state |j) = |5i/2,m = 1/2). 



|20] taking the form ([T]) by exciting both ions simultaneously with a bichromatic laser 
beam with frequencies uj± = ujq ± 6 where ujq is the qubit transition frequency and S 
is close to a vibrational mode of the two- ion crystal with frequency u (see figure [T]( a)). 
Changing into an interaction picture and performing a rotating-wave approximation, 
the time-dependent Hamiltonian 

H{t) = hn{e-''' + e^5*)e^'?(--'+'^^'=-')(a« + ) + h.c. (2) 

is well approximated by 

H{t) = -hrjnia^e'^"-^'^' + ae-'^''~^^')Sy (3) 

in the Lamb-Dicke regime where the Lamb-Dicke factor r] satisfies the condition 77 -C 1. 
In ([3]), we use a collective spin operator Sy = a'^^ + and denote the laser detuning 
from the motional sidebands by z/ — 5 = e. The Rabi frequency on the carrier transition 
is denoted fi, and a, are the phonon annihilation and creation operators, respectively. 
This Hamiltonian can be exactly integrated [21] yielding the propagator 

U{t) = D{a{t)Sy) exp {i{Xt - x sm(et))5j) , (4) 

where a{t) = ^(e*^* - 1), A = r]^n^/e, x = r/^^^Ve^ and D{a) = e""^""*'^ is a 
displacement operator. For a gate time tgate = 27r/|e|, the displacement operator 
vanishes so that the propagator ?7(tgate) = exp(iAtgate5'^) can be regarded as being 
the action of an effective Hamiltonian 

H,s = -hXSl = -2hX{l + ay® ay) (5) 

taking on the form given in ([T]). Setting Vt = |e|/(4?7), a gate is realized capable of 
maximally entangling ions irrespective of their motional state. In the description of 
the gate mechanism given so far, a coupling of the light field to the carrier transition was 
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Figure 2. Effect of non-rcsonant excitation of the carrier transition, (a) For C ~ Oj 
the gate starts at a maximum of the intensity-modulated beam. In this case, a Bloch 
vector initiaUy centered at the south pole of the Bloch sphere performs oscillations 
that are symmetric around the initial position, (b) For ( = 7r/2, the gate starts at 
the minimum of the intensity modulation. In this case, the average orientation of the 
Bloch vector is tilted with respect to its initial position. 

neglected based on the assumption that the Rabi frequency Q was small compared to 
the detuning 6 of the laser frequency components from the transition. In this case, small 
non-resonant Rabi oscillations that appear on top of the gate dynamics are the main 
effect of coupling to the carrier transition. Since a maximally entangling gate requires 
a Rabi frequency Q oc ri'H'J^^^, the question of whether Q 6 holds becomes crucial 
in the limit of fast gate operations and small Lamb-Dicke factors. Our experiments [9] 
are exactly operating in this regime, and it turns out that non-resonant excitation of 
the carrier transition has further effects beyond inducing non-resonant oscillations |12j . 
This becomes apparent by interpreting terms in the Hamiltonian in a different way: 
The red- and blue-detuned frequency components E± = £'ocos((co'o ± S)t ± () of equal 
intensity may be viewed as a single laser beam E{t) = E^+E^ = 2Eq cosiuot) coa^dt+Q 
that is resonant with the qubit transition but amplitude-modulated with frequency 5. 
Here, the phase which determines whether the gate operation starts in a maximum 
(C = 0) or a minimum {C, = vr/2) of the intensity of the amplitude- modulated beam 
has a crucial influence on the gate. This can be intuitively understood by considering 
the initial action the gate exerts on an input state in the Bloch sphere picture shown 
in figure [2J For short times, coupling to the sidebands can be neglected which justifies 
the use of a single-ion picture. The dynamics is essentially the one of two uncoupled 
qubits. The fast dynamics of the gate is induced by excitation of the ions on the carrier 
transition. For C = 0, the Bloch vector of an ion initially in state ||) will oscillate with 
frequency 6 along a line centered on the south pole of the Bloch sphere. For ( = 7r/2, 
the oscillation frequency is the same, however, the time-averaged position of the Bloch 



5 



vector is tilted by an angle 

= sm C (6) 

with respect to the initial state ||). This effect has a profound influence on the gate 
action. A careful analysis of the gate mechanism [12] taking into account the non- 
resonant oscillations reveals that the Hamiltonian Q is changed into 

Hit) = -nriVLi.a^e'^"-^^' + ae-'^''-^^')Sy^^ , (7) 

where 

Sy,ip = Sy COS ip + S;^smip , (8) 

and that the propagator (jlj) needs to be replaced by 

U{t) = exp{-zF{t)S,)D{a{t)Sy,^,) exp (z(At - x sin(et))52^) , (9) 

where the term containing F(t) = ^ {sm{6t+C) —sin () describes non-resonant excitation 
of the carrier transition. The dependence of the propagator on the exact value of ( is 
inconvenient from an experimental point of view. To realize the desired gate, precise 
control over ( is required. In addition, the gate duration must be controlled to better 
than a fraction of the mode oscillation period because of the non-resonant oscillation. 
Fortunately, both of these problems can be overcome by shaping the overall intensity 
of the laser pulse such that the Rabi frequency Q{t) smoothly rises within a few cycles 
2-11/6 to its maximum value legate ~ k 1/(4^7) and smoothly falls off to zero at the end 
of the gate. In this case, the non-resonant oscillations vanish and ([6]) shows that the 
operator Sy^^{t) adiabatically follows the laser intensity so that it starts and ends as 
the desired operator Sy irrespective of the phase (. For intensity-shaped pulses, the 
propagator (jl]) provides therefore an adequate description of the gate action. 



2.2. AC- Stark shifts 

In the description of the gate mechanism given so far the ion was treated as an ideal two- 
level system. AC-Stark shifts are completely insignificant provided that the intensities 
of the blue- and the red-detuned frequency components are the same since in this case 
light shifts of the carrier transition caused by the blue-detuned part are exactly canceled 
by light shifts of the red-detuned light field. Similarly, light shifts of the blue-detuned 
frequency component non-resonantly exciting the upper motional sideband are perfectly 
canceled by light shifts of the red-detuned frequency component coupling to the lower 
motional sideband. 

For an experimental implementation with calcium ions, we need to consider 
numerous energy levels (see figure [I](b)). Here, the laser field inducing the gate action 
causes AC-Stark shifts on the qubit transition frequency due to non-resonant excitation 
of far-detuned dipole transitions and also of other 6*1/2 ^ -D5/2 Zeeman transitions. The 
main contributions arise from couplings between the qubit states and the 4p— states 
that are mediated by the dipole transitions 51/2 ^ -P1/2, 'S'1/2 ^ -P3/2, -D5/2 ^ -P3/2- 
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Other transitions hardly matter as can be checked by comparing the experimental 
results obtained in [22] with numerical results based on the transition strengths [23] 
of the dipole transitions coupling to the 4p— states. For suitably choosen k-vector and 
polarization of the bichromatic laser beam, these shifts are considerably smaller than 
the strength A of the gate interaction. 

AC-Stark shifts can be compensated for by a suitable detuning of the gate laser. 
An alternative strategy consists in introducing an additional AC-Stark shift of opposite 
sign that is also caused by the gate laser beam [22j. This approach has the advantage of 
making the AC-Stark compensation independent of the gate laser intensity. In contrast 
to previous gates relying on this technique [21] where the AC-Stark shift was caused by 
the quadrupole transition and compensated by coupling to dipole transitions, here, the 
AC-Stark shift is due to dipole transitions and needs to be compensated by coupling to 
the quadrupole transition. 

For ions prepared in the ground state of motion {n = 0), a convenient way of 
accomplishing this task is to perform the gate operation with slightly imbalanced 
intensities of the blue- and the red-detuned laser frequency components. Setting the 
Rabi frequency of the blue-detuned component to fib = ^(1 + ^^e one of the 
red-detuned to fir = fi(l — ^), an additional light shift caused by coupling to the carrier 
transition is induced that amounts to ^i^^ = 2(fi^ — fi^)/*^ = — 8fi^^/5. Now, the 
beam imbalance parameter ^ needs to be set such that the additional light shift exactly 
cancels the phase shift (p = ^ac^gate induced by the dipole transitions during the action 
of the gate. Taking into account that tgatc = 27r/e and fi = |e|/(4?7), this requires 
e = (5Wkl)(0/vr). 

Apart from introducing light shifts, setting C, ^ also slightly changes the gate 
Hamiltonian ([5]) from ifcff = — AS*^ to H^s = —^{Sy + ^^5*^) [25] since the coupling 
between the states ||J,), |tT) is proportional to 2fibfir = 2fi^(l— ,^^) whereas the coupling 
between ||t), |||) is proportional to fi^ + fi^ = 2fi^(l + ^^). However, as long as ^ <^ 1 
holds - which is the case in the experiments described in the next section - this effect 
is extremely small as the additional term is only quadratic in C,- 

Another side effect of setting ,^ 7^ is the occurrence of an additional term oc Sza^a 
in the Hamiltonian. It is caused by AC-Stark shifts arising from coupling to the 
upper and lower motional sideband which no longer completely cancel each other. The 
resulting shift of the qubit transition frequency depends on the vibrational quantum 
number n and is given by 5^^^^ = {d>'rfVl?'/e)^n = {e/2)^n. Simulations of the gate 
action based on ([2]) including an additional term oc Sz accounting for AC-Stark shifts 
of the dipole transitions and power-imbalanced beams show that the unwanted term 
oc Sztt'^a has no severe effects for ions prepared in the motional ground state as long as 
^ <^ 1. However, for ions in Fock states with n > 0, this is not the case. Taking the 
parameter set = 0.05, u = (27r) 1230 kHz, and e = (27r) 20 kHz as an example, the 
following results are obtained: applying the gate to ions prepared in | li)\n = 0), a Bell 
state is created with fidelity 0.9996. For n = 1, the fidelity drops to 0.992, and for n = 2 
to even 0.976. This loss of fidelity can be only partially recovered by shifting the laser 
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Figure 3. Trap geometry and optical qubit liglit field generation. AOM 1 controls 
the overall frequency uji and the amplitude of the laser beams at 729 nm. AOM 2 is 
used to switch beam 2 which illuminates both ions simultaneously. When supplied 
with two frequencies Wrf ± S where a;rf/(27r) = 80 MHz, it creates a bichromatic light 
field in the first order of diffraction. AOM 3 switches beam 1 which addresses only 
one of the ions. Both beams are guided with single mode optical fibers to the trap. 
The geometric alignment of polarizations, trap axis and magnetic field is as sketched. 
Radio-frequency signals supplying the AOMs are indicated as dotted lines. 



frequency by S^^^\ the resulting fidelity being 0.996 and 0.988, respectively. For higher 
motional states, the effect is even more severe and shows that this kind of AC-Stark 
compensation is inappropriate when dealing with ions in a thermal state of motion with 
n ^ 1. Instead of compensating the AC-Stark shift by imbalancing the beam powers, 
in this case, the laser frequency needs to be adjusted accordingly. 

3. Experimental setup 

Two ^^Ca"*" ions are stored in a linear Paul trap with an axial trap frequency v/{2ti) = 
1.232 MHz corresponding to an inter-ion distance of 5 yum. The '^^Ca"'" optical qubit 
consists of the metastable state |t) = \D^/2,i^ = 3/2) with a lifetime of 1.17s and 
the ground state ||) = \Si/2,m = 1/2) (see figure [D^b)). These two energy levels are 
connected via a quadrupole transition at a wavelength of 729 nm. Laser light at 397 nm 
is used for Doppler-cooling and state detection on the 5*1/2 ^ P1/2 transition with an 
additional repumping laser at 866 nm on the -D3/2 ^ P1/2 transition. Fluorescence light 
is detected by means of a photomultiplier tube. For two ions we discriminate between 
ITT); {|ti) oi' lit)} and III), the populations of which are labeled by po, pi, and p2 
according to the number of ions fluorescing. 

A titanium sapphire laser [26], whose frequency is stabilized to a high finesse 
Fabry-Perot cavity [27], is used to excite the quadrupole transition for sideband- 
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cooling, frequency-resolved optical pumping and performing quantum logic operations. 
Frequency drifts of maximally 3Hz/s induced by the reference cavity are canceled by 
an automated measurement routine referencing the laser frequency to the optical qubit 
transition frequency and detecting the magnitude of the magnetic field at the ions' 
location of about 4 Gauss. 

The setup for controlling the laser driving the qubit transition is depicted in figure [31 
Laser light of 729 nm is sent to the ions from either of two directions, each beam having a 
maximum light power of 50 mW. Only when single-ion addressing is required we use laser 
beam 1 focused to a beam waist of 3 /im at the trap center with a /c-vector perpendicular 
to the trap axis and a polarization that couples to all possible transitions. All other 
operations are accomplished with laser beam 2 whose /c- vector encloses a 45 ° angle with 
the axis of the ion string and is perpendicular to the quantization axis defined by the 
direction of the magnetic field. With a beam waist of 14 /j,m at the trap center, this beam 
is adjusted to illuminate both ions with equal intensity. The polarization of this beam is 
set such that the coupling is maximal for Am = ±1 transitions whereas it vanishes for 
all other transitions. The amplitude, frequency and phase of both beams is controlled by 
the acousto-optical modulator (AOM) 1 which is driven by a versatile frequency source 
(VFS). Amplitude pulse shaping is achieved with a variable gain amplifier controlled 
by a field programmable gate array. All radio-frequency sources are phase-locked to an 
ultra-stable quartz oscillator. By triggering each experimental cycle to the AC-power 
line we largely reduce distortions caused by the ImG magnetic field fluctuations at 
50 Hz. 

The AOMs 2 and 3 are used as switches for laser beams 1 and 2 applied from 
different directions. The bichromatic light field with frequencies u!± = uq ± 5/(27r) 
is created by driving AOM 2 simultaneously with the two frequencies curf ± 6, where 
6 = u — e. A frequency difference of 25/(27r) = 2.4 MHz leads to a diffraction into 
slightly different directions with an angular separation as small as 0.025° such that the 
coupling efficiency to the single mode fiber is reduced by about 15% compared to a single 
frequency beam where AOM 2 is driven with Uj-f. To generate the collective 7r/2-pulses 
needed for analyzing the gate action, AOM 2 is driven with a single frequency Uj-f. A 
more detailed description of the apparatus is given in [4]. 

4. Measurement results 

The coupling strength of the laser to the qubit is calibrated by recording resonant Rabi 
oscillations on the qubit transition. In case of short gate operations, the large intensities 
lead to big AC-Stark shifts and saturation of the gate coupling strength [12] which in 
turn necessitate a fine-adjustment of the laser frequency and power. 
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Figure 4. Effect of amplitude pulse shaping on non-resonant population transfer 
caused by a bichromatic light field non-resonantly exciting the carrier transition. 
Experimental results are presented for a gate duration of tgate = 25 /xs. A comparison 
of the evolution of the populations po{k),pi{o),p2{») for a square-shaped pulse (a) 
with an amplitude-shaped pulse (b) shows a suppression of the strong non-resonant 
oscillations for the latter case. The slopes are shaped as a Blackman window with a 
duration of 2.5 ^s, the figure inset showing the definitions of pulse and slope duration. 
Numerical simulations suggest that the actual pulse shape is not so important as long 
as the switching occurs sufSciently slowly. The solid lines are calculated from ^ 
and (gl). 



4-1. Amplitude pulse shaping 

The merits of amplitude pulse shaping were studied by observing the time evolution of 
the populations pi at the beginning of the gate operation when the population transfer 
is dominated by fast non-resonant coupling to the carrier. Figure H] (a) shows the 
population evolution for the first 5 /is of a 25 /is gate operation based on a rectangular 
pulse shape. Averaging over a randomly varying phase (, we observe strong oscillations 
with a period of 27r/5 = 0.84 /xs. Panel (b) shows that the non-resonant excitations 
vanish completely after application of amplitude pulse shaping with a slope duration of 
2.5 /is corresponding to three vibrational periods of the center-of-mass mode. The slopes 
were shaped as a Blackman window [28], where the form of the shape is chosen such 
that a shaped and a rectangular pulse of the same duration have the same pulse area 
(see inset of panel (b)). Different pulse lengths are achieved by varying the duration of 
the central time interval during which the laser power is constant. The solid lines in the 
figure are calculated from ([9]) and (jl]) and fit the data well. 

4-2. AC- Stark shift compensation 

The AC-Stark shift caused by bichromatic light with spectral components each having 
a Rabi frequency of Q/{27r) = 220 kHz (for tgate = 25 /is) is measured by scanning the 
global laser frequency using AOM 1. The resulting populations after a gate operation are 
depicted in figure [5] (a). We observe a drop of the population pi to zero at a detuning of 
(27r) 37 kHz from the carrier transition. At this setting the ions are maximally entangled. 
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Global frequency detuning (kHz) Global frequency detuning (kHz) 

Figure 5. (a) Populations po{k),pi{o),p2{») after a single gate operation (tgatc = 
25 fis) where the global frequency detuning of the bichromatic entangling pulse is 
varied by scanning AOM 1. A maximally entangled state is achieved for a global 
frequency detuning of (27r) 37 kHz relative to the qubit transition frequency due to 
AC-Stark shifts, (b) Introduction of a beam imbalance ^ = 0.08 shifts the pattern 
of the populations by the required amount to fully compensate for the AC-Stark-shift 
(note the different x-axis offsets in (a) and (b)). The solid lines are calculated by 
solving the Schrodinger equation for the Hamiltonian given in ([2]) amended by a term 
accounting for the measured AC-Stark shift. 
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Waiting time t„ (|js) 

Figure 6. Population evolution of |||)(») and |],J,)(a) when scanning the waiting 
time between two 25 /is gate pulses in a Ramsey-like experiment. For this scan the 
detuning e was set to (27r)40kHz. In this set of data, the AC-Stark shift was only 
partially compensated by imbalancing the power of the two frequency components. 
From the sinusodial fits shown as solid lines, we infer an oscillation period of 258(4) /is 
corresponding to a residual AC-Stark shift of (27r)1.94(3) kHz. 

By changing the relative power of the bichromatic field's frequency components such 
that ^ = 0.08 the AC-Stark shift is compensated. This translates the whole excitation 
pattern in frequency space as can be seen in figure [5] (b). 

A more sensitive method to infer the remaining AC-Stark shift 6ac after a coarse 
pre-compensation consists in concatenating two gates separated by a waiting time 
in a pulse sequence akin to a Ramsey- type experiment [29] and scanning r^. This 
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procedure maps 5ac to a phase cj) = Sac^w which is converted into a population change 
P2 = cos^(</)), po = sin^(0) by the second gate pulse. For two ions, the corresponding 
Ramsey pattern displayed in figure [6] shows oscillations of the populations po and p2 
with a frequency of two times the remaining AC-Stark shift. 

4-3. Gate analysis 

A full characterization of the gate operation could be achieved by quantum process 
tomography [30]. At present, however, the errors introduced by single ion addressing 
and individual qubit detection are on the few percent level in our experimental setup 
which renders the detection of small errors difficult in the entangling operation. Instead, 
we characterize the quality of the gate operation by using it for creating different Bell 
states and determining their fidelities. 

For the Bell state = + 2|tT), the fidelity is given by F = (^i|p^^P|^i) = 
(p^l^P ^ + p^''^)/2 + Imp^'^i^, with the density matrix p"^^^ describing the experimentally 
produced state. To determine F, we need to measure the populations p2 +Po at the end 
of the gate operation as well as the off-diagonal matrix-element To determine the 

latter, we apply a 7r/2 pulse with optical phase to both ions and measure (ai^VF^) 
for the resulting state as a function of 0. This procedure is equivalent to measuring 
oscillations of the expectation value TT{P{(j))p'^^^) of the operator P{4>) = a^^^^a^^ where 
<^<j> = <^x COS (p + (Ty sin (p (see figure [9] (b) and (d)). The amplitude A of these oscillations 
equals 2|p^^|-|.| and is obtained by fitting them with the function Pm{<P) = A sin(20+0o)- 

Previous measurements [9] using |||) as input state have demonstrated Bell state 
fidelities as high as 0.993(1) (see figure [9] (a) and (b)) for gate times of 50 /is or 61 trap 
oscillation periods. By doubling the detuning to e/(27r) = 40 kHz we reduce the gate 
duration to only 31 trap oscillation periods and observe Bell states with a fidelity of 
0.971(2) which is remarkable considering the small Lamb-Dicke parameter of rj = 0.044. 
The detrimental effects illustrated in figure|4](a) are sufficiently suppressed by amplitude 
pulse shaping. 

Moreover, for a gate time of 50 /is the analysis was extended by applying the gate to 
the state |J,t) which is prepared by a 7r/2 rotation (beam 2) of both ions, followed by a 
TT phase shift pulse on a single ion performed with the far-detuned focused beam 1, and 
another tt/2 rotation applied to both ions. This pulse sequence realizes the mapping 

lii) ^ li + T)li + T) ^ li - T)li + T) ^ lit) (10) 

to the desired input state for the gate. Imperfections of single-ion addressing lead 
to an error in state preparation of 0.036(3). For the Bell state analysis, we measure 
the population pi to infer p^^u ~^ P^ll\v Unfortunately parity oscillations cannot be 
introduced by a collective pi/2 pulse acting on the state |t|) + «|it)- Instead, we 
transform this state into |tt) + ^lii) by repeating the steps of sequence ffTOl) as for the 
state preparation and measure again the coherence by performing parity oscillations. 

Figure [7] shows a comparison of the fidelity of the gate starting either in HI) or |||). 
The fidelity of a Bell state created by a single gate starting in |J,|) is 0.95(1). Taking 
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Figure 7. Bell state fidelities after gate operations applied to the input states m)(A) 
and |iT)(*) foi' ^gatc = 50^s. Taking into account the error for state preparation of the 
input states ||t) a-nd a similar error to measure the parity signal, we conclude that the 
gate operation works on all tested input states similarly well. The solid lines reflect a 
Gaussian decay of the parity fringe amplitudes as a function of the number of gates 
and a linear decay in the desired populations caused by the spectral impurity of the 
laser. For both input states the gate operation implies errors of less than 0.2 after 21 
consecutive applications. 

into account the errors for state preparation and the Bell state analysis we conclude that 
the entangling operation works equally well for Hf) as an input state. This hypothesis 
is supported by the observation that for both states we obtain a similar decay of Bell 
state fidelities with increasing gate number. 

Compared with our earlier results [5] where multiple gate operations were induced 
by varying the duration of a single bichromatic pulse, here we applied up to 21 individual 
amplitude- shaped pulses. Splitting up a long pulse into many shorter gate pulses has no 
detectable effect on the fidelity of the Bell states produced, and in both cases we obtain 
a Bell state fidelity larger than 0.80 after 21 gates. 

4.4- Gate errors 

As discussed in [9], the two dominant sources of gate errors are laser frequency noise 
and variations of the laser-ion coupling strength. Imperfections of the laser's frequency 
spectrum lead to incoherent carrier excitation and thus to a loss in the Bell state fidelity 
of 2 X 10~^ per gate. Coupling strength variations of 6Q/Q ~ 1-4(1) x 10~^ are the major 
cause for the Gaussian decay of the parity oscillation amplitudes. 

An error that was not investigated in detail before is the dependence of the Bell 
state fidelity on the global laser frequency detuning from the qubit transition frequency. 
Experimental results are shown in figure [HI The solid line fitting the data is calculated by 
numerically solving the full Schrodinger equation for different global frequency detunings 
and evaluating the fidelity. A second order frequency dependence of —9.6(3) x 10~^ Hz^ 
is found from calculations at the maximum point. This suggests that our laser's typical 
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Global frequency detuning (kHz) 

Figure 8. Fidelity as a function of the global frequency detuning of the bichromatic 
light pulse for from the carrier transition (here, the sideband detuning was set to 
e/(27r) — 20kHz). A maximum fidelity of 0.995(4) was found for a detuning of 
500 Hz from the transition center due to a residual AC-Stark shift. The solid line 
is obtained by numerically solving equation ([21) and taking into account the AC- 
Stark shift compensation by different powers of the blue and the red laser frequency 
component. At the maximum, the solid line experiences a second order frequency 
dependence of -9.6(3) x 10"^ Hz^. 

mean frequency deviation of 160 Hz contributes with 3 x 10~^ to the error budget. 

A further error source arises when the bichromatic beam couples to both ions with 
different strengths. By recording Rabi oscillations simultaneously on the two ions we 
conclude that both ions experience the same coupling strength ^2 to within 4%. From 
numerical calculations we infer an additional error in the measured Bell state fidelity of 
less then 1 x lO"'^. 

Another possible error source is heating of the COM-mode during the gate operation 
since the gate is not insensitive to motional heating in the parameter regime of our 
implementation. Using the calculation performed in [21], we find a fidelity reduction of 
AF = r/itgatc/2 where is the heating rate of the COM-mode. As in our experiments 
Vh = 3s^\ the fidelity is reduced by AF ^ lO"'^ for tgatc = 50 fis. 

5. A M0lmer-S0rensen gate with ions in a thermal state of motion 
5.1. Formal description of the time evolution 

In theory, the M0lmer-S0rensen gate does not require the ions to be cooled to the 
ground states of motion since its propagator (jll) is independent of the vibrational state 
for t = tgatc- For t ^ tgatc? however, the interaction entangles qubit states and vibrational 
states so that the qubits' final state becomes dependent on the initial vibrational state. 
Therefore, it is of interest to calculate expectation values of observables acting on the 
qubit state space after applying the propagator for an arbitrary time t. As will be 
shown below, simple expressions can be derived in the case of a thermally occupied 
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motional state. For the following calculation, it is convenient to define V{t) = exp{i'ySy) 
where 7 = At — xsin(et). We are interested in calculating the expectation value of the 
observable O given by 

0{t) =TT{OU{t)pM ® PAUit^) 

= Tr(pM ® pAD{-aSy)V^OVf){aSy)) 

= J2PnT^^AiPA{n\Di-aSy)V^OVDiaSy)\n)). (11) 

n 

Here, pM = ^nPn\^){^\ with pn = (w+l)" describes a thermal state with average 
phonon number n and Tr^ denotes the trace over the qubit state space with p^ being 
the initial state of the qubits. For two ions, the state-dependent displacement operator 
D{aSy) is given by 

D{aSy) = Po + P2D{2a) + P^2D{-2a), 

where Pa is the projector onto the space spanned by the eigenvectors of 5*^ having the 
eigenvalue A, with Pq = 1 — jSy, and P±2 = ^{Sy±2Sy). This decomposition allows for 
tracing over the vibrational states in (fTTj) since {n\D{a)\n) = exp(— |ap/2)£„(|Q;p), 
where denotes a Laguerre polynomial. For taking the trace, we note that 
^^p„(?T,|D(a)|?7,) is proportional to the generating function g[x,(3) of the Laguerre 
polynomial [21] given by 

n=0 ^ ^ 

Therefore, 

Y,Pn{n\b{a)\n) = -±-g{-^, \a\') exp(-|a|V2) = e-H'(^+^). (12) 
^-^ n + 1 n + 1 

n 

Using the abbreviation = VOV^ and (|T2l) . the expectation value 0(t) is given by 
0{t) = Tta{Ov{Ao + A4e-^l'^l'("+^) + Ai6e-^''l"l'("+^)}) (13) 

with 

Ao = PoPaPo + P2PAP2 + P^2PaP-2 
Ai = P2PAP0 + P0PAP2 + P~2PaPo + PoPaP-2 
Aie = P-2PAP2 + P2PAP-2 ■ 
For the initial state pa = \ |, one obtains 

A, = {Sl + Sl)/lQ, A4 = -5,/4, A,, = {Sl-Sl)/lQ 

and 

0{t) = 1^,a{Ov{{SI + ^ - 45,e-^l"l'(^+^) + [Si - 5,2)e~^^l"l'(^+^)}). 
lo 

To calculate the time evolution of the qubit state populations starting from state \[[) 
at t = 0, use of the relations 

e-'^'^yS.e'^^y = cos(47)^, - sin(47)^{^,, 
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yields the following expressions for the qubit state populations: 

p^(t) = -(3 + e-i6|ap{n+i) + 4cos(47)e-^l"l'("+5)) 
8 

= ^(l-e-^'^l"l'("+^)) (14) 

p,(t) = ^(3 + e-^'^l"l'("+5) -4cos(47)e-^l"l'("+5)) 
8 

The formalism presented here could also be used to calculate the contrast of a parity 
scan for thermal states of motion. In this case, the parity operator is given by 
P = 5*^/2 — 1. The 7r/2 carrier pulses transform this operator into an operator 
= (cos 05^ + sm(j)Syf/2 - 1. 

5.2. High fidelity Bell states of ions in a thermal state 

While many theoretical papers discussing M0lmer-S0rensen and conditional phase gates 
put much emphasis on the possibility of entangling ions irrespective of their motional 
state by using these gates, there has not been any experimental demonstration of this 
gate property up to now. The reason for this is that independence of the motional 
state, as predicted by (j4]), is achieved only deep within the Lamb-Dicke regime whereas 
experiments demonstrating entangling gates on hyperfine qubits usually have Lamb- 
Dicke factors on the order of r/ = 0.1 — 0.2 [101 El [32]. Therefore, all previous 
experimental gate realizations used laser cooling to prepare at least the motional mode 
mediating the gate in its ground state with n = 0. 

Figure [9] (a) illustrates the population evolution induced by the gate pulse for 
ground-state cooled ions initially prepared in the qubit states |J.|), figure [9] (b) displays 
parity oscillations for the produced Bell state. The corresponding time evolution and 
parity oscillations for ions that are merely Doppler-cooled to a thermal state with 
n = 20(2) are shown in figure [9] (c) and (d) respectively. As the coupling strengths 
on the upper and lower motional sidebands scale as oc \/n + 1 and oc ^Jn., non-resonant 
sideband excitation transfers population much faster from | [i,n) into | J,t,n ± 1), 
111,77, ± 1) as compared to the case of ions prepared in the ground state with fi = 0. 
After the gate time tgatc = 50 /xs, however, the undesired population pi nearly vanishes 
as in the case of ground-state cooled ions and the Bell state is again created. In 
the experiment, we find a population pi = 0.015(1) in the undesired energy eigenstates. 
The parity oscillations have an amplitude of 0.964(2), resulting in a Bell state fidelity of 
0.974(1). The reasons for the somewhat reduced fidelity as compared with ground-state 
cooled ions are currently not well understood. In part, the fidelity loss arises from a 
variation of the coupling strength on the vibrational sidebands as a function of n caused 
by higher-order terms in rj. However, for a thermal state with n = 20 and i] = 0.044 
calculations show this effect amounts only to additional errors of 7 x 10"'^. 

As mentioned in Section [2^2} the AC-Stark compensation by unbalancing the power 
of the red and blue frequency component is not applicable to ions in a thermal state. 
Instead, the laser frequency needs to be adjusted to account for AC-Stark shifts 6ac, a 




Figure 9. Measured population evolution for po(A)jPi(o),P2(*) and parity oscillations 
with (a.b) and without (c,d) ground state coohng. In the latter case, population is 
transferred faster into |Ti,n), |it,n) as compared to sideband cooled ions due to the 
higher coupling strength to the sidebands. In (c), the solid lines are a fit to the data 
points using (|14p with the mean phonon number n as a free parameter giving n = 20(2). 
The parity oscillations for the ions in a thermal state of motion have an amplitude of 
0.964(2). Combining both measurements, we determine the Bell state fidelity to be 
0.974(1). The data appearing in (a) and (b) are taken from [9]. Here the deviation 
of the solid lines from the data is caused by the AC-Stark shift compensation using 
C = 0.08. 

technique that works well as long as the AC-Stark shifts are smaller than the coupling 
strength A of the gate interaction appearing in ([5]) (otherwise, in the case 6ac ^ A, small 
laser power fluctuations give rise to large phase shifts). Therefore, care must be taken 
to choose the direction and polarization of the gate laser such that a favorable ratio 
\/Sac is obtained. In experiments with a gate duration of tgate = 50 /is, we achieved 
X/6ac ~ 3 and needed to shift the laser frequency by about 7.5 kHz for optimal Bell 
state fidelity. In future experiments, a further reduction of the AC-Stark shift could 
be obtained using a cr"'"-polarized laser beam incident on the ions along the direction of 
the magnetic field. In this geometry the AC-Stark shift is predominantly caused by the 
5*1/2 ^ -P3/2 dipole transition since the D^j2{rn = +3/2) state does not couple to any 
of the 4p Zeeman states. From calculations we infer a reduction of the shift to about 
2 kHz without compromising the gate speed. 

Fitting equations f|T^ to the population evolution data allows us to determine 
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the mean vibrational quantum number as n = 20(2). This value is consistent with 
independent measurements obtained by comparing the time evolution of the ions when 
exciting them on the carrier and on the blue motional sideband. 

6. Conclusions and outlook 

Until recently, entangling gates for optical qubits were exclusively of the Cirac-ZoUer 
type which require individual addressing of the ions. Compared to this type of gate the 
M0lmer-S0rensen gate gives an improvement in fidelity and speed of nearly an order of 
magnitude. The achieved fidelity sets a record for creating two-qubit entanglement 
on demand irrespective of the physical realization considered so far. Our results 
with concatenations of 21 of these operations bring the realization of more complex 
algorithms a step closer to reality. The implementation of a gate without the need for 
ground state cooling is of particular interest in view of quantum algorithms that require 
entangling gates conditioned on quantum state measurements that do not preserve the 
ions' motional quantum state. 

The optical qubit as used here is certainly not the best solution for long time storage 
of quantum information. Instead, qubits encoded in two hyperfine ground states whose 
frequency difference is insensitive to changes in magnetic field are preferable. These 
magnetic-field insensitive hyperfine qubits can store quantum information for times 
exceeding the duration of the gate operation presented here by more than four orders 
of magnitude [11 |3l H]. However, on such qubit states no high-fidelity universal gates 
have been demonstrated so far. Hence, our next experimental efforts will focus on 
implementing the M0lmer-S0rensen gate using ^^Ca^ ions. By mapping between the 
hyperfine qubit encoded in the ion's ground states and the optical qubit we will benefit 
from both of their advantages. 

Another interesting perspective of this gate is to create multi-qubit interactions 
between more than two qubits. A gate collectively interacting with all ions at the same 
time, in combination with a collective spin flips and a strongly focused off-resonant laser 
beam inducing phase shifts in individual ions, constitutes a basis set of Hamiltonians 
that offers the prospect of realizing complex multi-qubit operations such as a Toffoli-gate 
and a quantum error- correcting algorithm |33j . 
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